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We have formulated the problem of electron transport through interacting quantum dot system 
in the framework of self-consistent perturbation theory, and show that the current conservation 
condition is guaranteed due to the gauge invariant properties of the Green's functions and the 
generalized Ward identity. By using a generating functional for the statistics of the nonequilibrium 
system, we have obtained general formulae for calculating the current and the current fluctuations 
in the presence of arbitrary time-dependent potentials. As demonstration of application, we have 
studied the interaction effects on the finite frequency noise for electron resonant tunneling through 
an Anderson impurity, and obtained an analytical equation for the interaction effect on the finite 
frequency current noise within the Hartree approximation, which is an extension of the previous 
results obtained by Hershfield on zero frequency shot noise. 

PACS numbers: 73.23.-b, 73.63.Kv, 72.10.Bg 
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I. INTRODUCTION 

The transport properties of mesoscopic conductors 
have attracted wide research interest, and are of great 
importance for future nanoscale electronic applications. 
One usually probes the dynamics of electron in the con- 
ductors by measuring the dc and ac conductance in the 
linear response regime or driving the systems to the non- 
linear out-of equilibrium case. For instance, Gabelli et 
alii measured the ac conductance of a mesoscopic RC cir- 
cuit at GHz frequency, and found the violation of Kirch- 
hoff's law of impedance addition. Recent experimental 
advance will enable to probe the electronic processes in 
these systems in the high frequency region reaching the 
intrinsic time scales of the electron dynamics, therefore 
new interesting physical properties are expected to be 
observed i£ 

Some important progresses in this field are based on 
a wisdom that more information about electron dynam- 
ics can be obtained by measuring the current fluctua- 
tions or the higher current moments in these systems? 3 -* 4 . 
The theory of full counting statistics (FCS) of elec- 
tron current, which describes the probability distribu- 
tion of transmitted charge during a fixed time interval 
through a mesoscopic conductor, were developed within 
the scattering formulation^— and the Hamiltonian for- 
mulation for FCS were constructed based on the An- 
derson impurity model4~— It was also shown that the 
real part of ac conductance is related to the asymmetric 
parts of the frequency dependent current noise by a non- 
equilibrium fluctuation-dissipation theorem.! 11 : 12 There- 
fore, at present there are intensive research actives on the 
study of FCS problems both in quantum dot systems 1 ^— 
and diffusive conductors! 16 ' 17 

For the study of the time-dependent and finite fre- 
quency transport properties of mesoscopic systems, 
Blanter and Biittiker 4 emphasized the importance of con- 
sidering electron interaction effects, and pointed out that 



the current measured in frequency dependent transport 
experiment is a sum of the particle current and displace- 
ment current, while the displacement current is zero in 
the static case, but is essential for finite frequency trans- 
port. A theoretical approach that can address the fi- 
nite frequency electron transport and take into account 
the various electron correlation effects as well as the cur- 
rent conservation condition explicitly is highly worthy of 
invest igationi 18 ' 19 

A prototype mesoscopic system of nonequilibrium elec- 
tron transport with strong Coulomb interaction is a sin- 
gle quantum dot coupled to left and right leads, which 
can be described by the Anderson impurity model. The 
current through the quantum dofc2°- had been calculated 
by performing a second order perturbation theory of the 
Coulomb interaction strength U. However, the authors 20 
observed that the current is not conserved in this ap- 
proximation when the dot level is tuned away from the 
particle-hole symmetry point, and pointed out the ne- 
cessity of treating interaction terms by a current con- 
servation approximation. Hershfield 21 studied the cur- 
rent fluctuations in the Anderson impurity model sub- 
sequently, and obtained a formula for interaction effect 
on zero frequency shot noise in the Hartree approxima- 
tion. During the last decade, the shot noise and the 
FCS of current of the Anderson impurity model in the 
low-temperature Kondo regime have attracted a great 
deal of research efforts and interests, but we will re- 
strict our consideration only to the resonant tunnelling 
regime in this paper. For a noninteracting resonant tun- 
neling model, exact analytical formula of the finite fre- 
quency noise spectrum had been obtained^ The non- 
symmetrized noise spectrum of the resonant tunneling 
model had been studied recently within the scattering 
formulatio n 23 ' 24 , where some step and dip structures at 
finite frequency were shown, and the Hartree- Fock theory 
was applied to multi-level quantum dot system^ 

In the present work, we have formulated the theory of 
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electron transport through quantum dot based on the 
nonequilibrium self-consistent perturbation method^ 
which can guarantee the gauge invariance and current 
conservation condition for arbitrary time-dependent ex- 
ternal potential. We have studied the effect of the po- 
tential fluctuations in the quantum dot on the finite fre- 
quency noise within the Hartree approximation, and ob- 
tained an analytical equation for the interaction effect on 
the current noise, which is an extension of Hershfield's 
result to the finite frequency case. We will show that 
various correlation functions of nonequilibrium system 
including the vertex correction terms can be calculated 
in a systematic way by using functional derivations de- 
fined on the closed time contour and the external count- 
ing field method. Actually, the functional approach has 
been applied to a wide range of problems in mesosocopic 
system in the literatures, e.g. current and noise char- 
acteristics for single-electron transistors in the Coulomb 
blockade regim o 27 ' 28 or the Kondo regime^ photoas- 
sisted shot noise in mesoscopic conductor ^ charge trans- 
port to chaotic cavity 3 ^ and the time evolution of non- 
equilibrium quantum dot system^ etc. 

This paper is organized as follows: In section II, we dis- 
cuss the current conservation condition and the general- 
ized Ward identity for quantum dot based on the Ander- 
son impurity model. In section III, the formula for finite 
frequency noise spectrum including the interaction cor- 
rection term in the Hartree approximation is obtained. 
Section IV is devoted to the numerical calculations. In 
section V we summarize the results of this work. 



II. SELF-CONSISTENT PERTURBATION 
THEORY AND THE GENERALIZED WARD 
IDENTITY 

We consider the electrons transport through a single 
level quantum dot in the presence of external ac fields, 
the system will be described by the following Anderson 
impurity mode l 32 ' 33 

H = ^ £k n (t)cl va c kvi7 + ^ £d(t)dld a + U n d1 -n di 
+ ^2[t v e i ^(%l (T d a+ H.c] , (1) 

kr\a 

where 77 = L,R denotes the left and right leads, e^it) — 
f-kri + v v (t) and e d (t) — e d + v (t), with v v (t) and v (t) 
being the ac potentials in leads and in the dot, respec- 
tively. A^(i) is the gauge potential coupled to the tun- 
neling current from the lead rj to the dot. The strong 
Coulomb interaction term in the Hamiltonian prevent 
this model from an exact solution. But we can ap- 
proach this problem by doing perturbation theory on the 
Schwinger-Keldysh closed time path contour. Using the 
nonequilibrium Dyson equation, we write the equation 
of motion for the nonequilibrium Green's function of the 



quantum dot as follows 

[i^- t -e d (t)]G d(T (t,t') = 8{t,t') + J di 1 S(i,t 1 )G dCT (t 1) i / ) ■ 

(2) 

It should be emphasized that the time variables t and 
t' can be either on the forward or the backward branch 
of the time contour, and the integration over ti is also 
defined on the closed time path contour. The self-energy 
E(i, t') is obtained in the framework of self-consistent 
perturbation theory, and can be divided into two terms 

Z(t,t') = rt°\t,t') + Z u (t,t'), (3) 

where T,(°\t,t') is the dot level self-energy contributed 
from the tunneling between the leads and the quantum 
dot 

EW(t,0 = X)l t '»l a ^( t '^) e ~ i[ ^ W ~ A,(t,)+, ^ ,itl * , ' (tl)1 ' 

kij 

(4) 

with gkrj(t, t') being the bare Green's function of the lead 
without external ac potential field. Ej/(i, t') is the self- 
energy due to Coulomb interaction. In the self-consistent 
perturbation theory, it is a functional of the full Green's 
functions of the quantum dot. In order to illustrate the 
method of calculation, we consider only the first order 
approximation, and the interaction self-energy is given 
by the Hartree term 

X u (t,t') = U{n d *(t))5(t,t') . (5) 

We first study the gauge transformation properties 
and current conservation condition of the Green's func- 
tions of quantum dot. By making a transformation 32 

G dtT (t,t')_= G da (t,t r )e~ l ft dtlV0 ^\ the equation of mo- 
tion for G da (t,t') will be given by 

[i-^ - e d ]G dl7 (t,t') = 5(t,t') + J dtiE(t,ti)G ( fc r (ti,t'). 

(6) 

where the self-energy E(t, t x ) = T,^ ' (t,t') + 'E u (t,t') with 
E(°)(f,tO=El^| 2 ^(^Oe^ (t/) , (7) 

and in the Hartree approximation 

tu(t,t') = U{n ds (t))S(t,t') , (8) 

where the phase factor (j> v (t,t r ) = X tl (t) — X v (t') + 

L, dti[v v (ti) — vo(ti)], and G da {t, t') is a gauge transfor- 
mation invariant quantity. If one consider a gauge trans- 
formation: v (t) -t vo(t) + d t A(t), X v (t) -)• X n (t) + A(t). 
then it is easy to see that the phase factor <f) v (t,t r ), 
the self-energy E(t,t') and G da (t,t') are all gauge trans- 
formation invariant. Therefore, the Green's function 
G da {t,t') transforms as 

G da (t, t'; A) = e- lK ^G da {t,t')e^^ . (9) 
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The above gauge transformation is directly related to 
the current conservation condition in the quantum dot. 
Since under this gauge transformation, the change of the 
Hamiltonian to the first order of A is given by 



SH(t) = n d (t)d t A(t) + h(t)A(t) 



(10) 



where n^t) and j v (t) are the operators of the charge 
number in the dot and the tunneling current from the 
lead 77 to the dot, respectively. The gauge transformation 
invariance of the action leads to the continuity equation 



(11) 



In the out-of equilibrium steady state, the occupation 
number (rid(t)) is time independent, and the current con- 
servation condition X^nO'n) = is satisfied. 

Next we will follow the procedure as given in the Refill 
to give a derivation of the generalized Ward identity^ for 
this quantum dot system, which is closely related to the 
current conservation condition. In the spin degenerate 
case, the spin index a will be omitted. We consider the 
changes in the Green's function induced by the gauge 
transformation. From Eq.(9) the first-order change in G 



5G{t,t') = -i A(t)-A(t') G(t,t') , 
and it leads to the equation 



(12) 



5G < t - i '> a „A( ll) + E^A(' 1 ) 



Sv (h) 



SX v (ti) 



[A(t) - A(i') 



) G(t,t'), (13) 



where the functional derivatives of G can be denoted as 
Ao and A^, and they correspond to the time-ordered op- 
erator products as follows 

Ao(*,i';«i) = S -Pj^v = -(TcK(t)dUt')n d (ti)]) , 
dv (ti) 

(14) 

Wito - - -(2b[^(t)4(0^(*i)]> ■ 

(15) 

An integration by parts in Eq. (13), and demanding that 
the equation is satisfied for arbitrary A, straightforwardly 
leads to the well-known generalized Ward identity 

d tl A {t,t';ti) - £]A„(i,i';*i) 
v 

= i[S(t,h) - S(f \tx)]G{t,t') . (16) 

This identity leads to a relation between the vertex func- 
tions and the self-energy, which can be demonstrated ex- 
plicitly by introducing the following vertex functions 



r (M';«i) 



Sv Q (ti) 



(17) 



ft?- 1 (MO 
sx v (h) ' 



(18) 



One can see that these vertex functions are related to the 
time-ordered operators 



A a (t,t'-t 1 )= / <ft 2 disG(t,ta)r (t2,*3;*i)G(*3,*'). 



(19) 



A n (t,t';ti) = / dt 2 dhG{t,h)Tr,{t2,h)h)G{hJ) . 

(20) 

Thereby, the generalized Ward identity Eq. (16) can be 
rewritten in terms of the vertex functions 

d tl T (t,t';ti) - 5^r„(t,f;ti) 
n 

= i[«(f,t 1 )-5(t,ti)]G- 1 (t,* / ) .(21) 

This equation relates the vertex functions to the self- 
energy, since the inverse of Green's function G -1 = 
Gq 1 — £ is given explicitly as 

G-Ht,t') = [id t -e d -Vo(t)]5(t,t')-'£ ( - \t,t')-'E u (t,t') . 

The generalized Ward identity implies the gauge invari- 
ance and current conservation condition in this problem. 
Therefore it should be satisfied when we are investigat- 
ing the current fluctuations or time-dependent electron 
transport properties in this system by making approxi- 
mation calculation of the self-energy or the vertex func- 
tions. 



III. CURRENT FLUCTUATIONS AND THE 
HARTREE APPROXIMATION 

In this section, we study the current fluctuation and 
statistical problems in this quantum dot system. It is well 
known that the central quantity in FCS calculations^ is 
the cumulant generating function %(A) = J2q e iCJA P(Q), 
where A = (Ai, . . . , Ajv) are the counting fields and P(Q) 
is the probability for the charge Q = {Qi, ■ ■ ■ ,Qn) to 
be transferred through the respective channel^ during 
the measuring time T. For noninteracting electron sys- 
tem, the generating function is given by the Levitov- 
Lesovik formula^ within scattering matrix approach. It 
is observed that the cumulant generating function can be 
generalized to the system with time-dependent counting 
fields, and it can be written in terms of the nonequilib- 
rium Green's function defined on the closed time path 
contour as follows 

ln X (A) =TrlnG" 1 -TrE [/ G + $(G) , (23) 

where G _1 is the inverse of the full Green's function of the 
quantum dot given explicitly by Eq.(22). $ a functional 
potential constructed by summing over irreducible self- 
energy diagrams closed with an additional Green function 
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line^ The interaction self-energy Sj/ can be obtained 
from the functional 4> by 



Y,u{t,t') = 



5$ 



8G(f,t) 



(24) 



One can verify that when the counting fields A are as- 
sumed to be time-independent and the system is in the 
nonintcracting case, the above generating functional ar- 
rives at the Levitov-Lesovik formula. 

The electron current tunneling from the lead /; to the 
quantum dot can be obtained by a functional derivative 
of x(A) with respect to A^(i) 

it uw .e Slnx(X) 

{I " {t)) = l n-^W 

= -i e - J dhdhGihM)^ {t 2 M;t) • (25) 

Here the bare current vertex function (ta,ti;t) is 
given by 



5X v [t) 

= i[6{h,t) - 6(t 2 ,t)}^(t 2 ,h) .(26) 



Therefore the current through quantum dot is given 

„ „28,30 



(^(*)> 



dt x [G(t, i 1 )E(°) (t x , t) - E<°> (t, t x )G{t x , t)] . 

(27) 



By using the operational rules given by Langreth for con- 
tour integration^ it is not difficult to prove that this 
formula is exactly equivalent to the current formula ob- 
tained by Jauho et al. 37 for the time-dependent electron 
transport through an interacting quantum dot. 

We can introduce the interaction induced current 
vertex function which is related to the self-energy of 
Coulomb interaction 



(28) 



Then the vertex function defined in Eq.(18) is given by 



r„ (t 2 , h ; t) = ff (ta, h; t) + r£ (i 2 , *i ; t) 



(29) 



The current formula Eq. (25) indicates that the cur- 
rent is contributed solely from the bare current vertex. 
There is no contribution of interaction vertex correction 
to the current in this self-consistent perturbation ap- 
proach. However, we will show in the following that the 
interaction vertex indeed influences the current fluctua- 
tions. 

Next, we calculate the current-current correlation func- 
tions on the time contour and find that they can be rep- 
resented as the sum of two terms 



D m ,{t,t') = {T c 5I n {t)5Ir,-{t')) 



e 2 ^lnx(A) 
" h 5\ n (t)8Xn(t') 



to, 



(30) 



where the bare term is 



£©(*,*') = ^S w \G(t,t')^(t',t) + ^(t,t')G(t',t) 



dtidt 2 dt 3 dt 4 



G(h , fcjrg 5 (t a , t 3 ; t')G(t 3 , t 4 )r<°> (U, tv,t) 



and the interaction induced vertex correction term to the current correlation is given by 

.,2 



Dfy{t,t') = j I dhdt 2 dt 3 dt 4 G(ti,t2)r^(t2,t 3 ;t')G f (t3 ! t4)4 0) (t4,ti;t) 



(31) 



(32) 



Among the various current correlation functions, the correlation function for current noise is of particular interest, 
since the frequency dependent noise spectrum of current contains the intrinsic dynamics information of this quantum 
dot system. In a steady state without external time-dependent potential, the symmetrized noise spectrum (w) 
is given by the Fourier transform of the correlation function of current operators S vv >(t,t') = (SI n (t)SI v ' (t')) + 
(Sir)/ (t')SIri(t)) . It is noted that the correlation function for current noise can be written as 



S w (t, = D> (t, + D< (t, t') = (t, + S% (t, t') 



(33) 



where S^, (t, t 1 ) and S'fl (t, t') are contributed from the bare term and the interaction induced vertex correction term, 
respectively. 



The bare term S^j, (t, t') is obtained straightforwardly by using Langreth's analytical continuation rules^ In the 
absence of external ac potential, we can transform it to the frequency space, and express it in terms of the Green's 
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functions of quantum dot explicitly a a 38 ' 39 

-2 



dtu i 
~2t7 



^S w ir v [n v (uj 1 )G > (u>i +uj) - (1 -n^(wx + w))G < (w :L )] 

^|n n (o;i)(l - n n '(wi + w))G r (wi)G r (u;i + w) + n v >(ui)(l - n n (uJi + w))G a (u;i)G a (wi + w) 
+ [n^(wi)G Q (wi) - n^(wi)G r (wi)]G > (wi + uj) + G < (w 1 )[(l - n n >(cvi + ui))G r (u> 1 + u) 
-(l~n v (oj 1 +uj))G a (oj 1 +u)} - G < (uj 1 )G > (uj 1 + w)H + |w -> -wj . 



(34) 



In order to obtain the interaction effect on the noise spectra, we have to calculate the vertex function by functional 
derivation of the interaction self-energy with respect to the counting field: T^(ti,t 2 ;t) — -^xjjj) > where the inter- 
action self energy E[/(ti, £2) is given by Eq.(5) in the Hartree approximation. The technical details of our calculation 
is presented in Appendix B. After calculating the vertex function and transform it to the frequency space, we can 
obtain the interaction correction to the finite frequency current correlation function S^l (uj) as follows 



U 



UXnn 

u 



(0) 



(0) 



U 



\-U X r n ( n\uj) 



Or, 1 

u 



1-U X 



nv 



o,(0)/ \ 



o,(0)/ s. 



(35) 



where the various correlation and response functions are 
given in Appendix A. This equation is the central result 
of our paper. It is a generalization of the zero frequency 
noise result obtained by Hershfield 21 to the finite fre- 
quency case, and can be interpreted as the current noise 
contributed from the coupling of density fluctuations in 
the quantum dot to the current fluctuations. The first 
term in the Eq.(35) indicates that the correlation be- 
tween the density and the current S^2 (uj) can propa- 

gate forward in time via Xnn\w) and Xj'^n ( w ) to pro- 
duce current fluctuations. The second term represents 
the correlation Sj^ n (ui) propagates backward in time on 
the backward branch to produce current fluctuations. In 
the last term the correlation between the densities at the 
quantum dot Snn(&) propagates forward and backward 
on the closed time contour simultaneously and gives rise 
to current fluctuations. 



IV. NUMERICAL RESULTS 

To get better understanding of the interaction effect 
on the current noise spectrum for this quantum dot sys- 
tem, we will present some numerical calculations of the 
current noise at zero temperature. In our calculation, 
we take the coupling strength between the leads and 
quantum dot T as the units of the energy, and take the 
parameters ed = —1.0, U = 4.0 and the bandwidth 
D = 100. The applied bias voltage A/i = 3.0, with 
fXjj = — fin = Afi/2. For the system with symmetric cou- 
pling strength Tl = = T, we plot the current noise 
spectra for left and right leads in Fig. 1(a) and (b), re- 



spectively. Fig.l (a) shows that the bare noise spectrum 
5^(w) is positive and is an increasing function of the 
frequency. The interaction correction term S^ L (u>) can 
have negative values at an intermediate finite frequencies 
region. It is observed that the interaction correction for 
shot noise at zero frequency has a rather small positive 
value, which agrees with the previous result^ but the 
interaction correction becomes negative and more signif- 
icant when the frequency increases, and it goes to positive 
value again in the large frequency region. The maximum 
influence of interaction correction is obtained at the fre- 
quency which is largely determined by the energy differ- 
ence between the renormalized dot level Id — £d + U (rids) 
and the Fermi level of the leads. The sum of S^(uj) and 

SlI(u)) gives the noise spectrum after interaction correc- 
tion, which is a monotonously increasing function of the 
frequency. Fig.l (b) shows the noise spectra for the right 
lead (the drain side of this system). These noise spectra 
have more prominent features than that of the left lead. 
One can find a significant dip for the total noise spec- 
trum Srr(uj) at the frequency equal to the applied bias 
voltage (uj = A^=3.0). 

The various bare correlation and response functions 
utilized in the calculation of interaction effect of noise 
spectra are plotted in Fig. 2. As shown in Fig. 2 (a) that 

the density fluctuation spectrum Snn (w) of the quantum 
dot always has real positive values at finite frequencies. 
The real part of density response function Xnn(<^) is neg- 
ative at low frequency, which indicates that the screen- 
ing effect of electron decreases the Coulomb interaction 
strength U on the quantum dot. The imaginary part of 
Xnl (w) remains negative for all frequencies due to the an- 
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FIG. 1: (Color online) The noise spectra for quantum dot 
system in the symmetric coupling case, (a) for left lead, we 
plot the bare noise spectrum S^\(uf) (dashed line), the in- 
teraction correction term S^ L {ui) (dash-dotted line), and the 
noise spectrum after correction Sll{oj) (solid line); (b) the 
same for the right lead. Here, we take the parameters in the 
Anderson impurity model as = — 1.0, U = 4.0 in units of 
the coupling strength V, and we assume Fl = r_R = 1.0 for 
the symmetric case. 



alytical properties of the density-density response func- 
tion. The correlation and response functions between the 
density operator on the dot and the current operator for 
the left and right lead are plotted in Fig. 2(b) and (c), 
respectively. 

Fig. 3 shows the current noise spectra for asymmetri- 
cally coupled quantum dot system. Since the coupling 
strength Fl ^> Tr, we find that the magnitude of the 
current fluctuations in the right lead plotted in Fig. 3 (b) 
is much less than that of the left lead in Fig. 3(a), because 
of the tunneling rate between the right lead and quantum 
dot is much less than that of the left lead. The interac- 
tion correction terms also have negative value regions at 
finite frequencies both for the left and right leads. The 
noise spectrum in the drain lead (right lead) shows an 
evident dip structure at the frequency equal to the bias 
voltage. One can expect this kind of prominent features 
of noise spectrum can be detected in experiments. 



FIG. 2: (Color online) The spectra of various correlation and 
response functions involved in the calculation of the interac- 
tion effect on noise spectra in the symmetric coupling case (in 
units of e 2 /r). (a) The density correlation function Sira(w), 
the real and imaginary parts of the density response func- 
tion in the quantum dot. (b) The real and imaginary 
parts of the correlation function S^ n (uj) and response func- 
tion Xj'^'n ( LJ ) f° r t ne left lead, (c) The correlation function 
and response function of the right lead. The parameters used 
in the calculation are the same as in Fig.l. 
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V. CONCLUSIONS 

In this work, we have investigated the problem of elec- 
tron transport through quantum dot in the framework of 
nonequilibrium self-consistent perturbation theory and 
examined the current conservation condition. Based on 
the Anderson impurity model, we gives the current and 
current fluctuations formulae, which are valid in the pres- 
ence of arbitrary external time-dependent potentials, by 



FIG. 3: (Color online) The noise spectra for quantum dot 
system in the asymmetric coupling case . (a) For the left lead, 
we plot the bare noise spectrum S^l(io) (dashed line), the in- 
teraction correction term S^I(lo) (dash-dotted line), and the 
noise spectrum after correction Sll(u) (solid line); (b) The 
same for the right lead. We take the parameters in the An- 
derson model as = —1.0, U = 4.0 in units of the coupling 
strength F, and assume Fl = 1.0, Fr = 0.2 for the asymmet- 
ric case. 
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using nonequilibrium generating functional and the func- 
tional derivation method. We have calculated the in- 
teraction effect on the finite frequency noise spectrum 
of Anderson impurity model by taking into account the 
interaction vertex correction term with the Hartree ap- 
proximation, and obtained an analytical equation for the 
noise correction term at finite frequency, which corre- 
sponds to a generalization of the previous result on zero 
frequency shot noise. We have focused our attention on 
the symmetrized noise spectrum, one can expect that 
the nonsymmetrized noise spectrum and the ac conduc- 
tance can also be studied within the formulation pre- 
sented in this paper. We believe that the self-consistent 
perturbation theory on the Schwinger-Keldysh contour 
can lead to a unified approach to many interesting prob- 
lems in nonequilibrium electron transport through meso- 
scopic systems, and will give us deep understanding of the 
current fluctuation and energy dissipation phenomenal 
The functional method provides a convenient way to 
study the statistics of current fluctuations. We note that 
the Hartree approximation can only account for some 
interaction effects on the electron current noise in the 
resonant tunnelling regime. It is expected that future re- 
search work can treat the interacting effect beyond the 
Hartree approximation, and give us more information 
about the interaction effect on the out-of equilibrium dy- 
namics of electrons in the Coulomb blockade regime as 
well as the low-temperature Kondo regime. 
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Appendix A 

The response and correlation functions of density op- 
erator are given by 

XL(M') = -iO{t,t'){[n da (i),n da (t')]) , (Al) 

and 

S nn {t,t') = ({8n da (t),Sn da (t')}) , (A2) 

respectively. The current-density response function and 
correlation function 

X = -i0(t,t')(lj nt7 (t),n dr7 (t')}) , (A3) 

and 

S jvn (t,t') = ({Sj va (t) : Sn da (t')}) . (A4) 
By neglecting the vertex correction term induced by in- 
teraction effect, we can write the above functions in the 
frequency space explicitly as follows 



xSP-M - / £h)P-MG<<«. + «> + 0<MC {U1 + u „ , (A5) 

Sf,'H= f'!^\G<(u, l )G>{u, 1 +u) + G>(u 1 )G<(u, 1 +u,)\. (A6) 

/I 
{*T,n,(wi)[G r (« 1 )G r (wi + to) - G a {ui)G a {u)i - w)} + iT V G< ( Wl )[G r "(wi + w) + G a {Lo x - w)]} , 

(A7) 

S£l(w) = y^r j) {(l-n,(c 1 ))G''(c 1 )G < ( Wl + W )-n l; ( Wl )G I -( Wl )G > ( Wl +cG) 

+ n v (ui + uj)G > (uj 1 )G a {uj 1 + u) - (1 - n v (u)i + u))G < {cj 1 )G a (uj 1 + u) 

- G < (u 1 )G > (w 1 +uj)-G > (uj 1 )G < (u 1 + uj)} . (A8) 

In the Hartree approximation, the retarded/advanced Green's function G r / a (oJi) = l/[wi — e d — U{n dg ) ±i J2 V I\j/2]. 
The other Green's functions G < {u 1 ) = G r {w l )[iY ir ,T v n rt {w 1 )]G a (wi), and G > {u l ) = G r (o;i)[-^ r) r i( (l - 
n„(wi))]G a (wi). 
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Appendix B 

The interaction vertex function T^(t\,t 2 ;t) = ^sx^t)^ ' anc ^ tne se lf- ener gy in the Hartree approximation is given 

by 

3V(ti,t 2 ) = C^(n«»(ti)>5(ti,t2) = -iUGeituttWtu-h) . (Bl) 
In the spin degenerate case, we will omit the spin index cr. By using the following identity 

dt 2 dt 3 G(t 1 ,t 2 )[rW(i 2 ,t 3 ,t) + r^(t 2 ,t3,t)]G(t 3 ,t+) , (B2) 



■/■ 

we will derive the integral equation satisfied by the interaction vertex function. If we rewrite the vertex function 
F^(ti,t 2 ;t) = ^ {ti,t)5(ti,t 2 ), then the integral equation for T^(ti,t) is given by 

rV(t lt t) + iU J dt 2 G(t u t 2 )G(t 2 ,t+)fU(t 2 ,t) = -iUAW(t u t) , (B3) 

where 

4 0) (ti,t) - / dt 2 dt 3 G(h,t 2 )T^\t 2 ,t 3 ,t)G(t 3 ,tt) • (B4) 
Introduce a function M(t 2 ,t) which satisfies the following equation 

dt 2 [5(t u t 2 ) + iUG(t 1 ,t 2 )G(t 2 ,t+)]M{t 2 ,t) = 5(h,t) , (B5) 
then we obtain the solution for the vertex function as 

fjjf(ti,t) = -iU f dt2M(ti,t 2 )A^Ht2,t) . (B6) 

It is noticed that A^ ' (t 2 , t) is equal to the bare correlation function of the density operator on the quantum dot 

(<2,<) = {Tc5j n {t)5nds{t 2 )) without the interaction vertex correction as defined in Appendix A. The function 
M(t\,t 2 ) can be obtained by solving Eq.(B5). In the frequency space, it can be written in the matrix form of 

M+-(«) M++(c) J ~ I [/xl-' (0) M -1 + C/Xn„-' (0) (o;) J ' (B?) 

with A(cj) = [1 — Xn'n^ (<*>)] [1 — Xnri°'(^)]- Therefore, an explicit formula Eq.(B6) for the vertex function f^(ti,i) is 
obtained, which can be substituted to Eq.(32) to get the current noise interaction correction term. 
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